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Abstract 

Production cross sections for ti and ttj events at hadron colliders are cal- 
culated, including finite width effects and off resonance contributions for the 
entire decay chain, t —>■ bW —^ biu, for both top quarks. Resulting background 
rates to Higgs search at the CERN LHC are updated for inclusive H — > WW 
studies and for H ^ tt and H -^ WW decays in weak boson fusion events. 
Finite width effects are large, increasing ti{j) rates by 20% or more, after 
typical cuts which are employed for top-background rejection. 

I. INTRODUCTION 

tt production |[I[| is a copious source of VT-pairs and, hence, of isolated leptons at the Tevatron and 
the LHC. Top quark production will be intensely studied as a signal at these colliders. In addition, it 
constitutes an important background for many new particle searches. Examples include the leptonic 
signals for cascade decays of supersymmetric particles or searches for H — > W~^W~ [^|-Q and 
H^TT |,|,0] decays. 

Usually, ti production is considered in the narrow- width approximation (NWA), which effec- 
tively decouples top production and decay (see Fig. 0(a)). Whenever resonant top production 



dominates, this approximation is well motivated and it greatly simplifies matrix element evaluation 
and phase space generation. The NWA is also useful for single-resonant top production as shown 
in Fig. |l|(b) [|11|-|T^. In some cases calculations have been further simplified by also treating the 
decaying W bosons as on-shell particles. 
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FIG. 1. Feynman diagrams contributing to gg —>■ bbe'Uefi'^i^^g in the narrow-width approxima- 
tion: (a) on-shell ti production and (b) Wt single top production. The double bars indicate heavy 
quark propagators which may be treated as on-shell particles in various approximations. 



Naturally, the accuracy of these approximations needs to be tested, which requires a full cal- 
culation of off-shell effects. Restrictive selection cuts, as used for efficient suppression of ti back- 
grounds [^||, tend to be optimized against on-shell top production which may substantially enhance 
the relative importance of off-shell contributions. In applying the NWA, another problem inadver- 
tently arises: the Feynman diagrams in Figs. |l](a) and |l|(b) have identical initial and final states. 
When approximate cross sections, each specialized to a particular phase space region, are added 
up to obtain the total rate, double-counting can occur, and interference effects in overlap regions 
may not be properly accounted for. One thus needs a calculation which includes both resonant 
and non-resonant contributions, using finite width top-quark propagators, which correctly includes 
interference effects between the various contributions. The purpose of this paper is to present such 
a calculation for ti and tij production. In addition to merging resonant and non-resonant effects 



for the top quarks, we also include finite width effects for the W bosons, i.e. we consider general 
bbe~UefJ'^'^^i{j) final states. Finite top width effects, at the level of t ^ bW decays, have been consid- 
ered previously for tWb production |T^ . The complete set of Feynman graphs for pp —^ bW~^bW~ 
processes has been generated as well |T^ . However, the full treatment of lepton final states with spin 
correlations and off-shell contributions is new even for the ti case. To our knowledge, no previous 
calculation of finite width effects in tij production exists. 

The paper is organized as follows. In Section |T| we discuss various methods of including finite 
width effects and discuss their advantages and their practicality. We adopt the "overall factor 



scheme" and apply it in Section |T| to bbe Ug^^u^ production and the analogous process with one 



additional colored parton in the final state. While the matrix elements can, in principle, be generated 



with automated programs like MADGRAPH [|T^], the proper inclusion of finite widths, preservation 



of electroweak and strong gauge invariance, avoidance of double counting and of divergences in 



extraneous phase space regions requires manual intervention. In Section |T| we describe the content 
of our calculation, the various consistency tests, and other important features of the program which 
we have developed. 

Our program has already been used to study backgrounds to H —^ WW decays at the LHC [0. 
We expand on this analysis in Section ^ and use H -^ WW and H -^ tt decays, more precisely 
the backgrounds produced by ti and tij production, to exemplify the size of off-shell and on-shell 
contributions and compare our full simulations with previous background estimates. A summary 
and conclusions are given in Section 0. 

II. FINITE- WIDTH EFFECTS AND GAUGE INVARIANCE 

The inclusion of finite width effects is needed in order to avoid the singular behavior of the tree 
level propagators on mass-shell, p"^ — m? = 0. An approach that is straightforward to implement 
and, hence, well suited for automatic code generators like MADGRAPH/HELAS |[T^|T^ is to use a 



Breit-Wigner-type propagator with fixed width for all top and W propagators, making no distinction 
between time-like and space-like momenta. For massive fermions like the top quark one simply 
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FIG. 2. Dyson resummation of the imaginary part of the Wb contribution to the top self-energy. 
The resummed top propagator is represented by the cross-hatched blob. 

substitutes 

p — m p'^ — m^ p^ — m^ + tmi 

However, this technique, also referred to as fixed-width scheme, will lead to gauge-dependent 
matrix elements. Gauge invariance requires that the Ward identity 

k,v^ = -tg [tr\pi) - ^r\p2)] (2) 

(with pi = k + P2) be satisfied, where V^ = —ig'y^ is the gauge boson - fermion vertex. For the 
top-propagator of Eq. ([l|), the inverse, ^~ , even diverges at p"^ = m^ and the Ward identity is 
violated. The naive use of Breit-Wigner propagators does not produce consistent matrix elements. 
Calculating vertex functions and inverse propagators perturbatively, the Ward identity of Eq. (^ 
will be satisfied order by order. This is the basis of the so called fermion loop scheme for the W- 
propagator where, for a LO calculation, the imaginary part of the fermionic 1-loop corrections is 



included in both the vertex function and in the inverse propagator [1^^191 . For the propagator this 
corresponds to the Dyson resummation of the imaginary part of the W vacuum polarization. 

A theory driven solution of the finite width problem for the top quark propagator would gen- 
eralize this scheme. More specifically, a Dyson resummation of the imaginary parts of the 1-loop 
self-energy of the top quark, due to bW intermediate states (see Fig. H), results in the effective 
propagator 
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* ■p — mt + i^{p'^)-pPL Oip^ — {mw + rribY) 

Here Pl is the left-chiral projector and 
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In order to satisfy the SU{?)) Ward identity of Eq. @ one also needs to calculate the imaginary 
part of the ttg vertex (see Fig. |). We have checked by explicit calculation that this SU{?>) Ward 
identity is indeed satisfied. However, the effective ttg vertex already is too complex to be displayed 
here.Q For our applications we would need to know the imaginary contributions to ttgg and ttggg 
vertices as well. 
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FIG. 3. Effective ttg vertex including the imaginary part of the Wb loop-correction to the 
tree-level vertex. The imaginary part is obtained by cutting the triangle graph in all possible ways 
corresponding to on-shell intermediate states. 



A further complication arises from the need of electroweak gauge invariance in our calculations. 
Consider the simple process (or subgraph) depicted in Fig. ^: elastic scattering of a 6-quark and 



^It should be noted that the simplifications that led to concise results for the effective WW'j and 
WW'j'j vertices (see Refs. [1^,20, 21|) can not be applied when evaluating the effective ttg vertex. 



a longitudinal W boson. The longitudinal polarization vectors of the two W^s scale like vs/mw, 
which leads to a rise of the two subamplitudes with the center of mass energy: 

s [1 + ^7(5)] — m^ 

My,z ~ ^ • (4) 

t — my 

As is well known for the crossed process e~^e~ -^ W~^W~, gauge invariance of the electroweak 
couplings leads to a cancellation of these leading terms and results in partial wave amplitudes 
which do not grow with energy and, thus, respect partial wave unitarity p2|. When using the 



finite width top-quark propagator of Eq. (|]), with 7(3) ~ s at high energy, the width correction 
dominates the propagator at very high s and leads to A^j ~ const at high energy, thus spoiling the 
gauge theory cancellations between ^At and J^^y^z'- the bWi -^ bWi scattering amplitude violates 
unitarity in the J = \ partial wave at sufficiently high energy. The likely solution to this problem 
lies in adding the imaginary parts of btW, btWg vertices etc. to the loop scheme prescription, a 
solution which clearly becomes too involved for practical applications. 





FIG. 4. The electroweak gauge cancellation in bW^ -^ bW^ is perturbed by the resummed 
one-loop approach. 

When considering cross sections which contain contributions from both resonant and non- 
resonant amplitudes, 

a~/dPS.^,|M,,_+M/„e..P, (5) 



we need an alternative which guarantees gauge invariance and unitarity while allowing the effective 
substitution of propagators by a Breit-Wigner form in the resonant contributions, Aiij^res- In our 
work, we adopt the overall factor scheme [^ which starts from the observation that the zero width 
amplitudes, as derived from unresummed Feynman rules, provide a gauge invariant expression with 
proper high energy behavior. Multiplying the full lowest order amplitude by an overall factor 
{p"^ — m^)/(p^ — m'^ + irriT), for each resonant propagator, preserves gauge invariance but effectively 
replaces the tree level propagators, which are divergent on mass-shell, by finite-width Breit-Wigner 
propagators, 

2 2 2 2 

p - m^ p^-m^ 

Mif (r=0) -5 ^r~r~- F ~ ■^if,res. (F^O) + Mif^nonres. (r=0) -^ 5—— 7^ lOJ 

In the overall factor scheme, the cross section is thus calculated in terms of this modified am- 
plitude. Note that the additional overall factor is close to 1 in phase space regions where the 
non-resonant amplitudes yield significant contributions. Close to resonance, on the other hand, the 
amplitude effectively reduces to the dominant resonant terms. As we will show in Section III, this 
approach provides an excellent interpolation between on and off resonance regions, with ambiguities 
never rising beyond the 1-2% level. 

A second practical solution starts from the observation that the Ward identity in Eq. (0) remains 
fulfilled if one changes constant terms in the inverse propagator. This suggests another method to 
restore gauge invariance, which has been dubbed the complex-mass scheme P3,|23|. Its finite- 



width amplitude is derived from the full lowest-order amplitude with zero-width propagators by 
substituting all W, Z and top quark masses according to 



m 



Vm^ — imT (7) 



This scheme has recently been used in a single top quark study for the Tevatron [g^ . An unphysical 
consequence of the universal substitution (^ is that space-like propagators receive imaginary con- 
tributions, or that the top-Higgs Yukawa coupling, ht = rrit/v, receives an imaginary part. However, 
such effects are suppressed by factors T/m <^ 1 and would presumably not be noticeable in a LO 
Monte Carlo program. 



In our finite width ti and tij Monte Carlo programs, we have opted for the overall factor 
scheme. We expect very similar numerical results in the complex-mass scheme, which we have not 
implemented, however. 

III. CROSS SECTIONS FOR bbWW{j) PRODUCTION 

Unstable particles occur in several places in top quark pair production processes: in the form 
of decaying top-quarks and also as the W bosons which arise in their decays. In order to include 
off-resonance contributions for both, we are led to consider the full tree-level matrix elements for 
bbe~UefJ'^'^fj. final states ( "tt production" ) at 0{ala^), and final states with one additional parton 
("ttj production") at O(a^a^). 

A. Matrix Elements 

For pp or pp collisions, matrix elements for the following subprocesses need to be evaluated (we 
neglect CKM mixing): 

gg -^ bbe^h'efJ'^^ii , QQ -^ bbe^UefJ'^^fi for ti (8) 

gg -^ bbe'v^^^v^g , g^q^ -^ 66e~z/e/i+z/^g' , qq — > 66e~z/e/i+z/^5' for tij . (9) 

Representative Feynman graphs are shown in Fig. ^. Double resonant contributions include 
gluon radiation off initial state partons, the top quarks, and final state 6-quarks (Figs. |](a) and 
(b)). An example for a single resonant graph is shown in Fig. ^(c), while (d) depicts one of the 
non- resonant graphs. Electroweak gauge invariance of the bb -^ W^W~ subgraphs in (c) and (d) 
requires inclusion of 7 and Z exchange contributions. One such contribution is shown in Fig. ^(e). 
Others include ly-emission off the final state lepton lines (see Fig. ^ and discussion below). Our 
code includes finite 6-quark masses (set to a default value of mi, = 5 GeV). This allows to integrate 
over the entire 6-quark phase space, including the g ^ bb splitting region. A finite 6-quark mass 
necessitates new contributions, however, namely Higgs exchange diagrams like the one depicted in 



Fig. ^(f). Our code includes all these contributions. We avoid goldstone boson exchange graphs by 
working in the unitary gauge for the electroweak sector. 
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FIG. 5. Feynman diagrams contributing to gg —>■ bbe i'efJ'~^^fig with off-shell intermediate states: 
double- resonant (a,b), single- resonant (c) and non- resonant (d,e,f) contributions. 



TABLE I. Number of Feynman graphs contributing to different subprocesses. The numbers 
in brackets reflect the reduction of subamphtudes illustrated in Fig. || (see text). The columns 
correspond to different initial states. 
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Our calculation assumes different lepton flavors in the decay of the two PF-bosons. However, the 
amplitudes for this mixed lepton flavor case can also be used to obtain approximate results for same 
flavor processes, specifically e^UeC^i^e and fi^u^^^u^j, final states. The double- and single-resonant 
contributions (with respect to top) are identical for the mixed and same flavor sample. However, 
the latter also features (7,2' —>■ i~^i^) o (Z ^ i^ei^e) graphs, non-resonant from the view-point of 
top-decay, which do not occur in the mixed flavor case. Away from the Z-boson mass-shell and 
small £'^£~ invariant mass, these contributions are small, and the complete ifi^ cross section (with 
^1,2 = e, /i) is obtained by multiplying the result presented below with a lepton-flavor factor of 4. 
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FIG. 6. Amplitude factorization for 7 and Z decays. Explicit summation of the sketched helicity 
amplitude fragments leads to a significant computational reduction, as shown in Table |. 



The number of subamphtudes, corresponding to individual Feynman graphs, is sizable for the 
processes of Eqs. (H,P) and is listed in Table |. Constructing the code that evaluates the matrix 



elements was simplified by using automatically generated output of MADGRAPH |15| as a starting 



10 



point. However, the code had to be modified by hand in several ways. First, every Feynman graph 
amphtude needs to be multiphed with the proper overall factors 

2 2 

depending on its resonance structure. For ti production the factors correspond to a resonant top 
quark with momentum pi, + p^+ + p^^ (for which we use the shorthand notation Pb^+u ) and/or a 
resonant anti-top quark with momentum pie-p^- For tij production two additional factors ft{pb^i+v^,j) 
and ft{Pbe-Pej) appear, corresponding to gluon emission off the final state b or b. In both cases two 
W factors {fw{Pe-ue) ^.nd fwiPfi+u^,)) and one Z factor fzipe-v^n+v^,) need to be explicitly multiplied 
into amplitudes that are non-resonant with respect to these propagators. MADGRAPH generates 
all top, W and Z propagators automatically as resonant propagators, i.e. with a finite fixed width. 
These propagators can be viewed as the result of multiplying the zero-width tree-level propagators 
with the overall factor, or, alternatively, as obtained by the substitution 

[p\a : {p" - ml)-^ ^ (p2 -ml + imJ^aY' , (11) 

which we denote by the symbol \p\a in the amplitudes given below. For the first few Feynman 
graphs of Fig. ^ these changes can be summarized as 

b) M{[pb^,+ u^tAPhe-u,j\tAPe-v^wAP^^+v^w)* ft{p-he-u,)ft{Pbi^+v^j)fz{Pe-v,^+v^,) , (13) 

C) A^([P6M+^M]t'[Pe-pJiy, [P^+i.Jv^)*/t(P5e-pJ/i(P6M+i'Mi)/t(Pbe-p,i)/z(Pe-PeM+z^M) . (14) 

d), f) M{[Pe~u,]w, [Pi^+^M) * ft{Pbi,+ y^)ft{Pbe~9^ft{Pbi,+ u^j)ft{p-be-9,j)fz{Pe-u,t,+ u^) , (15) 

e) {Z) A^([Pe-pJiy, [P/.+ vJh-, [Pe-Pe/.+ i.Jz) * ft{Pbf,+ uJft{p-be-DjMPbf,+ u^j)Mp-be-p,j) ■ (16) 

In the overall factor scheme, imaginary parts for t-channel top-quark propagators are absent. The 
top-width was eliminated by hand from the MADGRAPH output for these space-like propagators. 
The effect of this modification should be small, however, since |p^ — ■m'^\ ^ mT if p^ < 0, due to 
T/m <^ 1. The number of Feynman diagrams for tij production is formidable, partially due to 
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(a) (5) 



FIG. 7. QCD radiative correction fragments, corresponding to (a) g ^ bb and (b) g -^ bbg 
splitting. The shaded blob represents the rest of the Feynman graph. Care has to be taken to count 
these contributions only once when combining relevant backgrounds. 

repeating sequences of subgraphs, like the ones depicted in Fig. ^. These subgraphs were combined 
to effective 7/Z-currents. As shown in Table 0, this procedure reduces the number of sub-amplitudes, 
which need to be calculated, by a factor of two or more.[| 

Since our calculation includes full matrix elements to a high order in perturbation theory, care 
has to be taken to avoid overlap with other backgrounds and double-counting. Closer inspection of 
the matrix elements for the ti and tij modes yields two groups of Feynman diagrams that are candi- 
dates for double-counting. They are schematically depicted in Fig. 0. Consider gu — > bbW~^W~u, a 
subprocess of tij production, as an example. Graph (a) represents g ^ bb splitting for a final state 
gluon. It constitutes a QCD radiative correction to gu —>■ W^W~ug and should be counted in the 
QCD WWjj background. For small bb invariant masses, graph (a) is a contribution to the gluon 
fragmentation, which must be counted only once. Group (b) features g — >• bbg splitting, where the 
6-quarks and one gluon are external particles. When the on-shell gluon is in the initial state and 



^Code generators have some freedom in the composition of basic elements when constructing 
helicity amplitudes, as explained in Section 2.7 of Ref. [|I6|. In rare cases, the algorithm employed 
by MADGRAPH is incompatible with the factorization outlined in Fig. |^. In these cases, suitable 
amplitudes were composed by hand. 
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the b is collinear to it, this group represents an as correction to bu -^ uW^W~h production. When 
both the h and the h are colhnear to the initial gluon, we are considering an a^ QCD correction 
to gu -^ uW^W~ . It is inappropriate to include either group as non- resonant contributions to 
tij production. These splitting contributions can be eliminated either through suitable cuts or 
by explicitly excluding the relevant Feynman graphs. A numerical comparison shows that typical 



selection cuts, like for the Higgs search to be considered in Section [I^, are sufficient to render the 



splitting contributions negligible. Since the analyses in Refs. |^ and IT^ do not include QCD cor 



rections to general WWjj production, we chose to omit these graphs in our final calculations. This 
is permissible since the graphs of group (a) and (b) necessarily contain electroweak interactions of u 
and d quarks. In the context of top backgrounds, the primary focus is on electroweak interactions of 
bottom and top quarks. When CKM mixing is neglected, the gauge invariance of our amplitudes is 
preserved when setting all electroweak interactions of the first two quark generations to zero. This 
procedure eliminates the graphs of groups (a) and (b) and avoids double-counting. 

When integrating over collinear g ^ bb configurations for initial state gluons, the differential 
cross section receives enhancement factors of order \og{fi'l/ml). Since as{'mt)\og{m^/ml) ^ 0.85, 
one may wonder whether our perturbative leading-order calculations are still reliable, or whether a 
resummation of these collinear logarithms is required. The correct treatment of the collinear region 
would include gb scattering, convoluted with the 6-quark parton density. Then, a subtraction of the 
gluon splitting term would also be required to avoid double-counting [p] , p!^ . 

In the cases at hand, the net effect of the 6-quark PDF contribution is small, however. For 
inclusive pp — > bbWW, only 10.4 pb are added to 622 pb (see Table 1 in Ref. |1^) . For the SM 



Higgs searches via weak boson fusion, discussed in Section |^, the tagging criteria and selection 
cuts are chosen such that for ti production both 6-quarks are resolved, with px > 20 GeV. This 
avoids the collinear region. For ti+1 jet production, the 6-quark observed as a forward tagging 
jet is required to have pt > 20 GeV, while the other 6-quark has no lower transverse momentum 
threshold. However, these collinear regions contribute little to the cross section within typical cuts. 
For example, in the H ^ tt -^ e^jji^^j. search with forward jet tagging cuts [0, the phase space 
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region with px < 20 GeV for the untagged b- (or b-) quark contributes only 8% to the total cross 
section. 

The smallness of these coUinear effects is related to the fact that we are generating jets and b- 
quarks oipT > 20 GeV as explicit partons in our calculations. In order to avoid double counting, the 
factorization scale in the 6-quark PDF should then be chosen as fij = 20 GeV, which mitigates the 
role of the collinear logarithms, aslog(yUf/m^).0 We conclude that a special treatment of collinear 
effects is not required in our studies. We regularize the 6-quark collinear region by the finite 6-quark 
mass, which provides a simplified but adequate model for the 6-quark PDF. 

Differential cross sections for top production in the narrow-width approximation are independent 
of rriH- However, once all off-shell effects are included, a dependence on ttih is caused by Higgs 
propagators that appear in the non- resonant contributions of Fig. Kf). This dependence has a 
neghgible effect on the rates considered in this paper. 

B. Phase Space Generator for tij Production and Numerical Tests 

While a simple phase space generator proved sufficient for ti production, for tij production with 
significant selection cuts a composite phase space generator that interfaces optimized mappings for 
double-, single- and non-resonant phase space regions becomes a necessity. These three regions are 
defined with the help of the two variables 

'mbw+-mt miw--mt 
nt = and Ut = (17) 

The double-resonant region is then defined by \nt\, \nt\ < ric, the single-resonant region by 
{\nt\ < ric A \nt\ > ric) V (|ni| > ric A \nt\ < ric) and the non-resonant region by |ni|, \nt\ > n^ (see 
Fig. H). In our calculations the boundary parameter Uc is chosen to be 8 by default (but see tests 



^As in any LO multi-parton calculation, the dangerous large logarithms are of the form 
aslog{mf/p'^j). An improvement of our calculation would first of all require the resummation of 
these contributions. Initial state collinear logs from g ^ bb splitting are unimportant by comparison. 
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FIG. 8. Phase space decomposition for bbWW final state based on the variables 
nt = {rnhw+ — T^t)/^t and rii = {mi^r- — mt)/Tt. The calculation presented here employs sep- 
arate, hand-optimized phase space generators for double- (D), single- {St and Sf) and non-resonant 
(A^) phase space regions. By default the boundaries are chosen at \nt^i\ = 8. 

below). Since the runtime of the tij program is fairly long, the code currently evaluates only the 
5*^ region and multiplies the result by two to account for Si- This symmetry assumption can easily 
be removed in the source code should the need arise. 

To assure the correctness of the programs three tests were performed. First, the Lorentz- 
invariance of the modified matrix elements is highly sensitive to errors. A suitable test variable 
can be defined in the following way: One generates a phase space configuration and evaluates the 
matrix element squared, summed over all helicity combinations. Applying an arbitrary boost (with 



^ Because of the massive external 6-quarks, the matrix element for any particular helicity combi- 
nation is not boost-invariant. 
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7 < 5) to all external momenta and re-evaluating the matrix element, one finally computes 

S \-Mboosted\ 



Evaluating 50-10^ phase space configurations (within the forward jet tagging cuts of Ref. pOd), 
we found an average (maximal) test variable of 2-10^^° (7-10^®) for tt and 1-10^^'^ (2-10^^) for tij. 
The test variable is highly sensitive to errors in the matrix element: Omitting the overall factors for 
a single amplitude that corresponds to a single-resonant Feynman diagram, increases the average 
and maximal test variable by 5-6 orders of magnitude. Second, the same test variable was used to 
check the modifications that are necessary to implement the reduction in Fig. p. 

A more recent MADGRAPH version can be used (with minor modifications) to automatically 
generate the full matrix elements of Eqs. (P) and (|^) with 8 and 9 external particles, respectively. 
These matrix elements can then be compared numerically to the factorized matrix elements in 
our programs. All overall factors need to be set to 1 for this test and finite widths for space-like 
top propagators need to be restored. A single, wrong coupling constant, e.g. gzu rather than gzd, 
increases the average and maximal test values by 11-13 orders of magnitude. 

Finally, the phase space generation for tt has been tested by comparing with known cross sections 
for the Tevatron and the LHC The composite phase space generation for tij has been checked by 
moving the boundary between phase space regions, i.e. changing the value of ric. For different cut 
sets and n^ = 4, 8 and 16 top-widths, we obtained results consistent within statistical errors of less 
than 1%. For tij one can explicitly take the narrow top-width limit and compare with the results 
in Ref. [§], for example. The programs passed all tests. 

In order to achieve an accuracy of 1% in a practical time, the programs use an enhanced version]^ 
of the VEGAS algorithm p7| , p8| . It applies importance sampling also to the summation of physical 
helicity combinations and separately optimizes suitably chosen combinations of subprocesses/phase 
space regions. 



The code is available at [http://hepsource.org/dvegas/ 
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TABLE II. Comparison of inclusive cross sections with and without finite top width effects. All 
results are given in pb. The first four columns represent our results: full ti cross section including 
off shell effects, contribution within \miyw+ — iTit\ > ricTt or \mi^^~ — mt\ > UcTt with Uc = 10 and 
15, difference to the NWA result, and this difference as a fraction of the total. The Wt column 
uses results from Table 1 in Ref. [^ for pp — > bbWW and Table 1 in Ref. for pp — > bbiuiu. The 



results in the second row include leptonic r decays, but not b decays. In Ref. [0 different methods 
are used for Wt to reduce overlap with ti: *^ Inibw — rnt\ > lOFj cut ^^ \mbw — rnt\ > 15Tt cut. 



process 


full 0{ala^'^) 


\mi)W — mt\ > lO/lSFt 


full - NWA 


rel.contr. 


Wt + wn 


rel.contr. 


pp -^ bbWW 
pp -^ bblviv 


622 
39.2 


63.9/49.6 
4.0/3.1 


25 

1.4 


+4% 
+4% 


66.4*/51.6§ 
4.0 


+ 11.1/8.6% 
+10% 



C. Numerical Results 



First numerical results were obtained using CTEQ4L parton distribution functions as a default, 
with0 as{mz) = 0.132 p9|. The renormalization and factorization scales firj are fixed to the top 
mass, nit = 175 GeV. In contrast to the Tevatron, top production at the LHC is dominated by the 
gluon-gluon channel and hence noticeably affected by uncertainties of the gluon density. To assess 
the impact of effects related to PDF and asirnz) choice, we compared CTEQ4L and CTEQ5L 
parton distribution functions [^. Studying the basic process pp —>■ ti, one obtains 622 and 510 pb 
for CTEQ4L and as{mz) = 0.132 and 0.118, respectively. On the other hand, for CTEQ5L and 
as(^z) = 0.127 and 0.118 one gets 560 and 487 pb, respectively, suggesting an overall uncertainty 
of about 10 or 20% related to the choice of PDF and as(mz) related to it. With selection cuts, we 



found similar deviations. For example, in Table ^ below, switching from CTEQ4L to CTEQ5L 
leads to a 4-8% decrease in cross sections at the forward tagging cuts level. 

A first comparison of our full calculation with single top production cross sections in the literature 



,(5) 



This as{mz) is based on Ag^^, as determined in the PDF set fit. 
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is presented in Table |T[ The two rows correspond to bbWW and bbi'^ui'i? final states, i.e. in the 

second row off-shell- H^ corrections are included in our calculation. In our simulation, bbe^ufi^i? 

results were multiplied by a factor 5.48 to account for all lepton flavor combinations, including 

leptonic r^ decays. Off-shell contributions from W~t + W~^i final states were calculated in Ref. |13[ 

with several definitions of the off-shell region. These results agree with ours at the 4% level (see 

second column). The agreement improves to the 1-2% level once double counting effects in the 

non- resonant regions (A^-regions of Fig. |^) are taken into account. Compared to this W~t + W^i 

cross section, however, the actual enhancement of the ti cross section due to off-shell contributions, 

calculated as the difference of our full result minus the cross section in the NWA (see column three), 

is only about half as large. 

The problem can be traced to the fact that a Breit-Wigner distribution has long tails, more 

precisely 

rim^nrf 1 _ ^ ^ ^ _ / _ ^X 

J{m-nT)2 (g2 - m2)2 + (mr)2 27r 2mT V rnr J ^ ' 

i.e. for each of the two top-quark resonances, about 2% of the cross section is located outside an 

n = 15 top-widths window, for example. This 4% resonant contribution, which is double-counted 

when combining the rate in the \mi)W — mt\ > l^Tt region with the ti cross section in the narrow 

width approximation, accounts for the difference between our "full minus NWA" result and the 

estimate in terms of the W~t + W^i cross section. This means that off-shell ti calculations which 

directly take W~t + W^i cross sections as the off-shell rate, tend to produce a serious overestimate 

of the number of extra events. 

Our method of choice to include finite-width effects involves overall factors that are expected 

to be small in non-resonant regions as explained in Section ||. In order to test this expectation 

in a realistic context, we compared the differential cross section obtained with overall factors with 

various proxies in different phase space regions (see Fig. §). In the non-resonant region A^ the proxy 

is the tree-level matrix element with unmodified tree-level top quark propagators, because they are 

not singular and a good approximation to the full propagator in this region. In the single-resonant 

regions St and St the proxy is derived from the proxy used in region A^ by replacing the potentially 
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singular top (in St) or anti-top (in St) propagators with fixed-width propagators, i.e. 

i{;^ + nit) 
p^ — ml + irritTt 



(20) 



100 




60 80 100 

PTj.mm [GeV] 

FIG. 9. Integrated transverse momentum distribution for the additional massless parton (jet) in 
off-shell tij production: aijp^^) = S^^A^'^ I^Vtj) dprj- No selection cuts are applied. For p™™ = 40 
GeV one obtains attj ~ Ctt = 39.2 pb (dotted line, see Table y). 



In the double-resonant region D the proxy matrix element is a subset of all Feynman diagrams, 
namely all diagrams with at least one time-like top propagator. Since all these top propagators are 
potentially resonant here, they all feature the fixed-width form of Eq. (|20|). Each region is covered 
by 20-28 bins. Bin sizes are adjusted so that each value roughly has the same integration error. For 
each bin the relative deviation from the proxy is estimated by 



da 



bin 



fac 



da 



bin 



proxy 



da, 



bin 



proxy 



(21) 



We studied this measure for tt production without cuts as well as the forward tagging cuts of 



Refs. and 0] and found similar results. In the double- and single-resonant regions the overall 
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factor scheme deviates from the proxies by 2% or less, and the deviation drops to less than 1% in 
the non-resonant region. We therefore estimate the error associated with our finite-width scheme 
to be of (9(1%), which is comparable to the statistical error of our results, but completely negligible 
compared to missing higher order QCD corrections. These results suggest that our method to 
include finite-width effects provides reliable results, not only for fairly inclusive cross sections, 
where finite width effects are strongly suppressed, but also when complex selection cuts result in 
on- and off-shell contributions of similar size. 

Comparing inclusive cross sections, we find Bau = 39.2 pb in Table |l[ Including an additional 
jet we obtain Bauj = 80 pb, when a common pt > 20 GeV cut is imposed on the additional jet. 
This raises the question, whether the cut should typically be chosen higher, such that auj ^ ctu, 
i.e. atij is sub-dominant, as typically expected for a NLO QCD correction. The integrated ptj 
distribution of the non-6 parton is shown in Fig. ^. For p™™ = 40 GeV one obtains auj ~ o"tj. Real 
parton emission cross sections saturating the LO cross section at low px of the extra parton are a 
sign of copious multi-jet production in actual data [pl|] . 

The additional parton emission is dominated by initial state radiation. This can be inferred 
from the invariant mass distributions of the potential top-quark decay products, the bW^ and the 
hW^j systems shown in Fig. [1^. The mi,y/+ invariant mass distribution in Fig. |iy(a) contains 85% 
of the total cross section in the displayed 165-185 GeV window around the top resonance (±6.4rt). 



In contrast, the same m^vF+j invariant mass window (insert of Fig. |TO](b)) accounts for only 10% 
of the total cross section. Final state radiation is relatively unimportant in ttj production at the 
LHC. 

IV. APPLICATION TO SM HIGGS SEARCH AT THE LHC 

An important application of tij production as a background occurs in Higgs physics. Higgs 
mass limits have recently been pushed to itih > 114.1 GeV by the LEP experiments |^2[. As 



a result, the LHC search for H — > WW [Q-|^ and H ^ rr PJ9|JTI]|| decays in the intermediate 
mass range has gained even greater importance. For the H -^ WW — * I'^t^^j. decay mode, top- 
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6000 




185 




1000 



(a) (b) 

FIG. 10. Invariant mass distribution of potential top-quark decay products (a) bW~^ and (b) 
hW^i in tij production, including full finite width effects. No selection cuts are applied, except for a 
Ptj > 20 GeV cut on the final state massless parton. Cross sections correspond to one lepton-flavor 
combination only. 

quark decays constitute the largest reducible background. The impact of this background has been 
analyzed, in the narrow-width approximation, for Higgs masses around 170 GeV [§,^§ and most 
recently for a light Higgs boson with tuh ~ 115 GeV ||^. This background also plays a role in 
the H ^ TT ^ e^/i^^y decay mode, which was analyzed in the narrow-width approximation in 
Ref. [|l^. In this section we present updated results for these analyses obtained with our parton- 
level Monte Carlo program that includes off-shell top and W effects and takes into account the 
single-resonant and non-resonant contributions.^] 



^The analysis in Ref. already features off-shell tt and ttj background estimates calculated with 
the programs presented here. 
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TABLE III. Numbers of expected pp -^ bbiuiu background events in the ATLAS H — > WW* 
Higgs search for an integrated luminosity of 30 fb^^. Basic ATLAS TDR background suppression 
cuts are apphed [1,0. ATLAS detector effects are simulated as described in Section II. F in Ref. 0. 
The left columns show our results for the four W -^ ev, ^v lepton combinations. The right columns 
are the corresponding results from Table 2 in Ref. [^. For both categories, the relative change due 
to the inclusion of finite width effects is given as an enhancement factor in the third column. 



rriT window [GeV] 


top in NWA 


full off-shell 


rel. chg. 


ti 


ti+Wt 


rel. chg. 


120-150 


50 


93 


1.9 


41 


206 


5.0 


130-160 


50 


96 


1.9 


66 


204 


3.1 


140-170 


43 


83 


1.9 


51 


146 


2.9 


140-180 


49 


95 


1.9 


61 


164 


2.7 


150-190 


35 


68 


1.9 


56 


111 


2.0 



A. tt Backgrounds to Inclusive H -^ WW Searches 



The ti background calculations (without an additional jet) are most relevant for inclusive H -^ 
WW searches [|]-|§], where Higgs production is dominated by the gluon fusion process. In Tables fT| 
and |IV| we compare the relative contributions from off resonant top effects for two selections of 
H -^ WW -^ I'^C^^j. events in ATLAS, as described in Refs. [§],§[. The selection looks for two 
isolated, opposite charge leptons of pt > 20, 10 GeV within the pseudo-rapidity range |?7^| < 2.5 and 
of invariant mass min < 80 GeV. Events must have significant missing Et, $rp > 40 GeV. A small 
angle between the charged leptons favors H —>■ WW decays versus backgrounds. Finally, a veto 
on additional jets in the central region is very effective against the 6-quark jets in the top-quark 
backgrounds. The main difference between Tables |Tl| and |3 is the definition of this veto on jet 



activity with px > 15 GeV in the central region. It is imposed within \rij\ < 3.2 in Table |T| and 



within |?7j| < 2.4 in Table |V. 

In addition, the selection looks for events inside the Jacobian peak of the dilepton-^^ transverse 
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FIG. 11. Transverse mass distribution for the low- luminosity ATLAS selection cuts of Ref. 0, 
Section 19.2.6, and Ref. 0. The dotted curve shows the distribution obtained when treating the 
top quarks in the narrow-width approximation while the solid curve includes all off-shell effects. 



mass distribution, as indicated in the first column of the two tables. Here the transverse mass is 
defined as 



mT{U,$^) = J2p^4$j,[l - cosM{a,$T)] . 



(22) 



Fig. [Tl| shows this transverse mass distribution underlying the ATLAS TDR analysis. Off-shell 
contributions raise the normalization of the ti background by about a factor of 2, but have little 
effect on the shape of the background, for itit ^120 GeV. The ATLAS analysis attempts to take 
(single-resonant) off-shell effects into account via on-shell Wt calculations. As the comparison in 
Table |T| shows, our unified calculation indicates a lower background increase than the combined 
on-shell tt and Wt calculations would suggest. This observation is consistent with the small increase 
due to off-resonant effects observed in Table |l|. 
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TABLE IV. Same as Table |T|, but for the optimized background suppression cuts of Ref. 0. 



niT window [GeV] 


top in NWA 


full off-shell 


rel. chg. 


ti 


ti+Wt 


rel. chg. 


120-150 


364 


536 


1.5 


107 


346 


3.2 


130-160 


363 


544 


1.5 


122 


306 


2.5 


140-170 


310 


465 


1.5 


82 


215 


2.6 


140-180 


352 


528 


1.5 


107 


254 


2.4 


150-190 


251 


374 


1.5 


87 


168 


1.9 



For a precise comparison with the ATLAS simulation, the nature of the programs (parton-level 
MC with full matrix elements vs. event generator with parton showers/hadronization) is currently 
too different. In particular our program does not allow to simulate the effect of the central jet veto 
on extra gluon radiation in the event. Also, we did not include any detector efficiencies. We would 
expect an extra suppression, by perhaps a factor of 3 to 4, of the "full off-shell" results in Tables |T| 
and |V| due to these effects, i.e. the apparent agreement of the ti cross section with our NWA result 
in Table |TTl| is somewhat fortuitous. This suspicion is confirmed by the fairly large disagreement 
between our NWA or full off-shell results and the PYTHIA simulation in Table [TVl . 

At present the source of these discrepancies is not completely understood. The rapidity distri- 
bution of 6-quarks in our full matrix element calculation appears to be wider than in the PYTHIA 
generated events, potentially due to approximations in the top decay chain in PYTHIA. This makes 
a veto over a smaller pseudo-rapidity range less efficient. On the other hand we cannot simulate 
the effect of additional gluon radiation. A reanalysis of these effects, combining full matrix ele- 
ments with a parton shower Monte Carlo, including hadronization, is clearly warranted, given that 
top-quark production constitutes about half the background to the inclusive H — > WW search 0. 

While the overall background normalization requires further study, the ratio of "full off-shell" 
and NWA results is expected to be robust, i.e. it will be little affected by detector efficiencies and 
by higher order gluon radiation. These ratios are given in the third columns, marked "relative 



change", in Tables HI and IV. Our results clearly indicate that the increase in top-backgrounds 
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due to the inclusion of off-shell effects (the Wt contribution in the ATLAS analysis) is substantially 
smaller than previously thought. 

In a LO calculation, substantial uncertainties arise from ambiguities in the choice of factorization 
and renormalization scales. For the double resonant phase space configurations, fir = P-f = iTit or 
the transverse energy of the top quarks are well-motivated choices. When forcing the 6-quarks to 
low pt values by a central jet veto and simultaneously enhancing the off-shell phase space regions, 
a smaller renormalization and/or factorization scale may be more appropriate. The impact of a 



lower scale on the full off-shell results in Table |T| is demonstrated by the scale choice 



mt if \mhw+ - mt\ < ATt A \miy^r- - mt\ < ATt 
fir J = { (23) 

20 GeV if \mhw+ - nitl > ATt V \mm^- - md > ATt 

where the value of 20 GeV is motivated by the veto threshold for central jets. Resulting cross 

sections are about 65% higher than the rates obtained with mt as a universal scale. A NLO 

calculation would be needed to distinguish the virtue of either choice. At present, this variation 

indicates the uncertainty of our LO results. 

B. ti{j) Backgrounds in Weak Boson Fusion Studies 

In addition to the inclusive search for H -^ WW events, weak boson fusion (WBF) presents a 
very attractive search channel for H -^ WW and H ^ tt events and will likely play a crucial role 
in the measurement of Higgs boson couplings to fermions and gauge bosons |3^ . Top quark decays 



again form an important background in these searches and, due to the central jet veto proposed for 
the reduction of QCD backgrounds, off-shell effects might be important. The signal and backgrounds 
for H -^ WW -^ e^fi'^'prp and H ^ tt ^ e^fi'^'pj. in WBF were analyzed in Refs. |^ and |10|, 



respectively, with the top-quark backgrounds determined in the NWA, however. With our new 
programs we are able to update these background estimates, including off-shell contributions. 

Because of the two additional jets which are present in the WBF process qq -^ qqH, the dominant 
top quark background arises from tij events. Event selection for WBF requires two tagging jets, 
of Pt > 20 GeV, which are widely separated in pseudo-rapidity (jr/i — ri2\ > 4.2) and which have a 
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very large dijet invariant mass, Mjj > 650 GeV. By definition, the ti background has both 6-quarks 
identified as tagging jets while in the tij background exactly one 6 or 6 is taken as a tagging jet. 
The two 6-quarks rarely have a large enough dijet mass or are far enough separated to satisfy the 
tagging criteria, and this leaves tij events as the dominant background to WBF. This result holds 
in both the NWA and with inclusion of off-resonant effects as is evident in both Tables |V] and [V^, 
at all cut levels. 

The veto of central jets of pt > 20 GeV is effective against the extra 6-quark jet and off-shell 
contributions are only slightly enhanced after this cut (see line "6 veto" in the Tables). Overall, 
off-shell contributions are fairly modest, increasing the NWA results by about 20% (see the second 
last columns in the Tables). This means that our new complete calculation of off-shell effects in 
tij production confirms the conclusions about the observability oi H ^ tt and H -^ WW events 
reached in Refs. ^ and |10|. Tables and ^ display our updated results. Precise definitions of 



cuts are given in the earlier papers. A breakdown into subprocesses and phase space regions of the 
overall tij background of 351 fb to if ^ WW -^ e^fj.^'^j., after forward jet tagging cuts, is given 
in Table IVTTl 



V. SUMMARY 

Top-quarks are a very important source of lepton backgrounds to new physics searches at the 
LHC: the large production cross section typical for a strong interaction process combines with a 
sizable branching fraction into leptons which, due to the large mass of the W, often survive isolation 
cuts. Suppression techniques for top quark backgrounds, like a veto on central jets, which is very 
effective against the 6-quarks of the t -^ bW -^ Mv decay chain, enhance the relative importance 
of off-resonant effects and may exacerbate errors introduced by approximate modeling of matrix 
elements. Severe cuts select the tails of various distributions and these phase space regions may 
well differ from the ones for which the models were optimized originally. General purpose Monte 



Carlo programs like PYTHIA ^^ or Herwig |^ should thus be gauged against matrix element 
programs. 
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TABLE V. ti{j) background cross sections for H -^ WW -^ e^fi^'^j, for rriH = 160 GeV in pp 
collisions at ^/s = 14 TeV. Results are given for various levels of cuts and are labeled by equation 
numbers from Ref. |^. All cross sections are given in fb. Cuts and other calculational details are 
described in Ref. [^]. The integration error is 1% or better. The signal over background ratio is 
also shown. Cross sections not listed here are as in Table I in Ref. Istl. 





top 


in NWA 






full off-shell 




cuts 


tt 


ttj 


S/B 


tt 






tij 


S/B 


forward tagging (10)-(12) 


12.4 


308 


wl/65 


13.0 


+4.4% 


351 


+14% 


«l/67 


+ b veto (13) 




43.5 


1/5.1 






51.4 


+18% 


1/5.6 


+ Mjj, angular cuts (14)-(16) 


0.0551 


4.67 


1.1/1 


0.0761 


+38% 


5.42 


+16% 


1.0/1 


+ real r rejection (17) 


0.0527 


4.34 


1.7/1 


0.0737 


+40% 


5.09 


+17% 


1.5/1 


Psurv,20 {x 0. 29) + (18) 


0.0153 


1.26 


4.6/1 


0.0214 


+40% 


1.48 


+17% 


4.2/1 


+ tag ID efficiency {xO.74) 


0.0113 


0.932 


4.6/1 


0.0158 


+40% 


1.09 


+17% 


4.2/1 



TABLE VI. ti{j) background cross sections for H ^ tt ^ e'^/x^^y for M^ = 120 GeV in pp 
collisions at ^/s = 14 TeV. Results are given for various levels of cuts and are labeled by equation 
numbers from Ref. |T0[. All cross sections are given in fb. Cuts and other calculational details are 
described in Ref. [|lOl. The integration error is 1% or better. The signal over background ratio is 



also shown. Cross sections not listed here are as in Table I in Ref. |jTO|. 





top 


in NWA 






full off-shell 






cuts 


ti 


tij 


S/B 


tt 




ti 


/■ 


S/B 


forward tagging (7)-(10) 


13.5 


357 


1/1100 


15.9 


+17% 


436 


+22% 


1/1100 


+ b veto (11) 




50.1 


1/550 






63.6 


+27% 


1/550 


+ ^j, (12) 


11.1 


43.0 


1/74 


13.2 


+19% 


55.2 


+28% 


1/83 


+ M,, (13) 


0.593 


12.9 


1/32 


0.712 


+20% 


15.8 


+22% 


1/34 


+ non-r reject. (14, 15, 17) 


0.00303 


0.257 


1/5.8 


0.00365 


+20% 


0.293 


+14% 


1/5.8 
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TABLE VII. Distribution of the tij background to H —^ WW -^ e^/i^jf^^, with forward tagging 
cuts, among subprocesses (labeled by initial partons) and phase space regions (see Fig. ||). The cut 
between double resonant (D), single resonant {St + St), and non-resonant regions (A^) is set at 8 
top quark widths. Cross sections are given in fb. 





D 


St + Si 


N 


99 


171 


28.6 


1.1 


gq + qg 


128 


18.8 


0.64 


qq 


2.92 


0.35 


0.009 



In this paper we have presented results for two new programs which allow to model the ti —>■ 
bbWW —>■ hM^vd,'!) decay chain at tree level, including full angular correlations of all top and 
W decay products and with proper interpolation between double-resonant, single- resonant and 
non-resonant phase space regions. These full correlations are available for ti and ttj production. 
Electroweak and SU{'i) gauge invariance is maintained throughout, by employing the overall factor 
scheme for the Breit Wigner propagators of all unstable particles in the t and i decay chains. 

Comparing to earlier calculations of off-shell effects, via the gh -^ Wtb production cross sections, 
we find excellent agreement when avoiding the top-quark resonance for the Wb system. However, 
some earlier combinations of ti and Wtb cross sections have involved substantial double counting, 
leading to an overestimate of backgrounds in e.g. Higgs search analyses at the LHC. In addition, 
the full simulation oi V — A couplings in the t -^ Wb -^ iub decay chain is available with our 
programs and may have sizable effects in background estimates. 

When considering ti backgrounds to Higgs searches, off-shell effects are most pronounced in 
inclusive H —>■ WW analyses, where both 6-quark jets may be vetoed. In weak boson fusion studies 
the additional jets in the signal selection make the background suppression due to the jet veto less 
severe, which diminishes the overall importance of off-shell contributions. 

Our programs now allow detailed study of these off-shell effects in top pair backgrounds with 
zero or one additional parton in the final state. 
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